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Model flocks in a steady vortical flow 
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We modify the standard Vicsek model to clearly distinguish between intrinsic noise due to imperfect 
alignment between organisms, and extrinsic noise due to fluid motion. We then consider the effect 
of a steady vortical flow, the Taylor Green vortex, on the dynamics of the flock, for various flow 
speeds, with a fixed intrinsic particle speed. We pay particular attention to the morphology of the 
flow, and quantify its filamentarity. Strikingly, above a critical flow speed there is a pronounced 
increase in the filamentarity of the flock, when compared to the zero flow case. This is due to the 
fact that particles appear confined to areas of low vorticity; a familiar phenomena, commonly seen 
in the clustering of inertial particles in vortical flows. Hence the cooperative motion of the particles 
gives them an effective inertia, which is seen to have a profound effect on the morphology of the 
flock, in the presence of external fluid motion. Finally we investigate the angle between the flow 
and the particles direction of movement and find it follows a power law distribution. 

PACS numbers: 47.32.Ef,47.63.-b,87.18.-h 


I. INTODUCTION 

Flocks of birds, schools of fish or swarms of insects are all 
examples of self-organised aggregations, offering some of 
the most spectacular displays of motion in nature. The 
complexity of the animals moving together, seemingly as 
one, has intrigued observers for centuries. Indeed, the 
ability of a flock of birds to change direction apparently 
simultaneously and spontaneously led Edmund Selous, 
who published numerous books on natural history, to be¬ 
lieve birds to be telepathic [I|. 

To mathematically model animal aggregations it is typ¬ 
ical to use models based on self-propelled particles, such 
as the Vicsek model 0. Here the direction of motion is 
typically based on the idea that flocking animals tend to 
align themselves to others within a local “radius of in¬ 
teraction” or based on a number of nearest-neighbours, 
as observed in starling flocks Q. More complex mod¬ 
els with additional short range repulsion have also been 
investigated 0, and a recent study extended the Vicsek 
model to include decision making, needed for example by 
birds to make a synchronised landing Q. 

We note that this is not the only approach with a num¬ 
ber of interesting studies based on nonlinear partial dif¬ 
ferential equations In addition a recent study by 
Ballerini et al. Q questioned the validity of approaches 
based on the standard Vicsek model, arguing that the ra¬ 
dius of interaction should be based on a topological dis¬ 
tance, rather than the metric distance used in [ 3 . How¬ 
ever we believe that the standard Vicsek model Q is 
the ideal starting point to investigate the effect of ex¬ 
ternal fluid motions on flock dynamics and morphology, 
and note that the behaviour of marching locusts could 
be modelled using a similar approach Q. There are also 
interesting extensions to model human crowd stampedes 

i- 

One of the drawbacks of the standard Vicsek model is 
that all sources of noise in the system arise in a single 
term, which is both temporally and spatially uncorre¬ 


lated. It is perhaps more natural to distinguish between 
two sources of noise, which we would expect to arise in 
such systems. The first is intrinsic noise, due to the fact 
that animals will never perfectly align. Here we see no 
issue with approach used in 0 where uncorrelated noise 
seems an appropriate approximation. However animals 
move in fluid environments, and so will be subject to ex¬ 
ternal noise, due to the motion of the fluid. In almost 
all natural flows a combination of large lengthscales and 
small viscosities mean that the flow is turbulent 0. 
Hence the external forcing will be subject to complex 
spatio-temporal correlations which are not present in the 
standard Vicsek model. 

A recent study by Khurana & Ouellette m sought 
to address this issue by considering a slight modification 
to the Vicsek model, in an unsteady flow made up of a 
summation of random Fourier modes, with a prescribed 
Energy spectrum consistent with the Kolmogorov spec¬ 
trum of classical turbulence. The focus of [TJ was the 
stability of flocks, particularly in the limit of a small ra¬ 
dius of interaction, where they showed an external forcing 
with spatio-temporal correlations was far more efficient 
in destabilising the flock. Here our interest is in the effect 
of fluid motion on coherent flocks, with larger values for 
the radius of interaction. 

One may be tempted, therefore, to surmise that it is 
then necessary to consider model flocks in fully turbu¬ 
lent flows. However, implementing a full direct numeri¬ 
cal simulation of the Navier-Stokes equation, necessary to 
produce a turbulent velocity field is a formidable compu¬ 
tational challenge [l^ . Moreover interpretation of results 
is not necessarily straightforward or unambiguous, given 
the strongly nonlinear nature of fully developed turbu¬ 
lence. Even the ‘reduced’ KS model considered in [ij 
is beyond what we wish to consider here, and it fails to 
reproduce many of the important features of true turbu¬ 
lence, such as intense vortical structures [T^. There is 
a long tradition in fluid dynamics of using simple ‘toy’ 
models, to probe the underlying physics of a problem. 
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FIG. 1: Snapshots of the Vicsek model in the statistically steady-state (left 77 = 0.1, R = 0.2, right 77 = 0.1, R = 0.8). Arrows 
indicate the particles direction of motion. 


with a certain level of control. This is the approach we 
take here, turning to the Taylor-Green vortex 0 as an 
ideal toy flow. Indeed this flow has been applied to fields 
as diverse as swimming microorganisms |l5| . the settling 
of inertial particles [la, [l^ and Magnetohydrodynamics 
0 . 

The aim of this paper is to investigate the effect of 
a simple vortical flow on a Vicsek flock, paying partic¬ 
ular attention to the flock morphology. We note that 
the shapes of schools of fish and flocks of birds have re¬ 
ceived both theoretical 0123 and experimental atten¬ 
tion 0, [2lj . However less attention has focused on the 
role that external fluid motion could play in affecting 
flock morphology. 

The plan of the paper is as follows: In Section m we 
describe the standard Vicsek model. We also introduce 
statistical measures of the shape and structure of the 
flock, namely Ripley’s K function and Minkowski func¬ 
tionals. In Section m we extend the model to include 
a smooth external forcing due to a steady vortical flow. 
Detailed comparisons of the statistical measures of shape 
and structure are made and we also study the angle be¬ 
tween the particles and the flow. We close with a discus¬ 
sion of the results in Section m and highlight possible 
directions for future studies. 


II. THE VICSEK MODEL 

We begin by introducing the Vicsek model [3, a well- 
known analogue of the Ising model of ferromagnetism. 
We consider N particles in a two-dimensional square box 
(with sides of size L), with periodic boundary condi¬ 
tions. Whilst the systems we are interested in are three- 
dimensional, by working in a low dimensional space we 
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FIG. 2: (Golor online) The time averaged Besag’s function 
{L{r)), Eq. ([5]) plotted as a function of spatial scale r for 
simulations with 77 = 0 . 2 . L{r) > 0 indicates clustering of 
particles at a given scale, L{r) < 0 indicates segregation. 


gain a huge computational speedup which allows for a 
thorough examination of parameter space. A follow up 
study to investigate flocking in three-dimensional time 
dependent flows is planned in the near future. Each par¬ 
ticle has a position Xi(t) and an intrinsic, self-driven, 
velocity All particles are assumed to move with 

the same speed, V, and a particles intrinsic velocity is 
determined by 

Vj = (Vcos(0i), Vsin(6»i)), (1) 

where 9i determines the direction the particle moves in. 
Key to the model is that 9i is periodically (at each time 
increment) determined from the average of the particle’s 
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own direction, plus the directions of its neighbours within 
a critical radius, i?, such that 


Bi = atan2 | ^ ^ sin6»j , ^ ^ cosBj | +?7Ci- 

\ |xi-X3|<_R |xi-X3|<_R / 

( 2 ) 

The final term in Eq. m is used to model both intrin¬ 
sic and extrinsic noise in the system. Specifically Q is 
a uniformly distributed random variable on the interval 
[—1,1] and rj is the intensity of the noise. From these sim¬ 
ple rules Vicsek et al. 0 observed the emergence of col¬ 
lective flocking or swarming of the particles. The global 
order of the system can be characterised by computing, 
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( 3 ) 


which is a convenient measure to establish if the astern 
has reached a statistically steady state. We follow and 
model a system of fV = 500 self-propelled particles in a 
domain of size L = 27r with a speed V = 1. Particles are 
evolved according to an Euler scheme such that 


Xj(f -1- At) = Xj(t) -1- AtVi{t) 


where at each tilmestep Vi{t) is updated according to 
Eqns. o & @. Note Q is drawn randomly at each 
timestep and we take At = 0.1. For large value of rj co¬ 
operative motion is not observed, as the system is domi¬ 
nated by noise, hence we restrict our parameter values to 
r] € [0, 0.5], R £ [0,1]. In order to be as computationally 
efficient as possible all distance evaluations are computed 
using a k-d tree (^ . 

Simulations are evolved for 200 time units. Within a 
short period of time (~ 20 time units) '4’v,R saturates 
to a statistically steady value, which for the low noise 
intensity considered here is close to or equal to unity, 
indicating substantial global alignment of the particles. 
Snapshots of the particles for simulations with rj = 0.1, 
R = 0.2, 0.8 are displayed in Fig. [U the effect that R 
plays in determining the structure of the flock is clear. 

In order to study this effect more quantitatively we use 
Ripley’s K function, a statistical pattern analysis method 
used as a measure of spatial clustering. For N particles 
this is defined as 


k{r) 


^ sp /(|xj - Xjl < r) 
N ^ N -1 

^^3 


( 4 ) 


where I is the indicator function (1 if the Euclidean sep¬ 
aration is less than r, 0 otherwise), and r G [0,L]. Note 
that this can easily be adapted to account for the periodic 
boundary conditions used in this study. Particle cluster¬ 
ing results in K{r) increasing faster than if particles are 
distributed in a spatially random manner, that is, if they 
follow a Poisson distribution. Ripley’s K function for a 
Poisson-distributed data set takes the form K{r) = Tik. 
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FIG. 3: (Color online) An example of a flock formed from the 
standard Vicsek model. Red (filled) points indicate where 
the separation between that point and its nearest neighbonr 
is less than ^/2, where I = L/y/N is the expected separation 
between random points. Blue (open) points are those which 
do not satisfy this criteria. The black structure shows the 
a-shape [ 2 J constructed from the red points, which is used 
to quantify the morphology of the flock. 


For a linear scaling of Poisson-distributed data, it is com¬ 
mon to normalise K{r) and define Besag’s function as 

L{r) = kK{r)lTT — r. (5) 

The advantage of Besag’s function is that it has a simple 
interpretation, it is zero for randomly distributed points, 
takes positive values for particles clustered over the spa- 
cial scale r and is negative if the particles are dispersed 
over a given scale. To increase our statistics we apply 
temporal averaging, and use angled brackets () to de¬ 
note temporal averaging in the statistically steady state 
(where '4>r],R has saturated); this is done every 10 time- 
steps. Figure [2] shows the averaged Besag’s functions for 
five different values of interaction radius R (rj = 0.2). 
A noticeable migration of the peak of (L) to the right 
is observed, as well as a broadening of the values of r 
where (L) > 0. Both are consistent with clustering of 
particles over larger scales with increasing R. We want 
to emphasise that this is not a particularly novel result, 
however it is important to establish ‘baseline’ statistics 
of the model, before studying the effects of flocking in 
the presence of a structured external flow. 

It is also instructive to introduce a statistic which mea¬ 
sures the characteristic size of the flock. Whilst there is 
no unique way to do this, we find it convenient to mea¬ 
sure 

D = argmax(L(r)), (6) 

r 

the location of the maximum value of the time averaged 
Besag function. The variation of (D) will be investigated 
later when we introduce a modified Vicsek model to ac¬ 
count for an external flow. 
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FIG. 4: (Color online) The Taylor-Green vortex flow. The 
pseudocolour plot displays the flow’s vorticity, Eq. black 
arrows indicate the velocity held, Eq. ([5|). 


We also study the morphology of the flock, and we are 
not aware of any study which has previously done so. 
To do this one must first compute the shape of the flock. 
Again there is no unique way to do this and it is clear that 
results will be somewhat dependent on the method that 
is used. Here our goal is not to perform a detailed analy¬ 
sis of flock morphology, for example to compare with that 
observed in bird flocks, insect swarms etc. But instead 
to quantify the change in morphology when an external 
flow acts on the particles. Hence we choose a compu¬ 
tationally efficient method, and apply this consistently 
in all simulations. We use the particle positions to de¬ 
fine an a-shape [ 2 ^ (a generalisation of the convex hull). 
Our procedure is the following: Given N particles in a 
domain of area the expected separation between ran¬ 
dom points is £ = L/y/N. Hence we define a subset of 
points Yi = {xi : jx^ — x,/1 < jt}, where i' denotes the 
closest particle to x,;. We take 7 = 0.5, which leaves the 
set of points whose distance to the nearest neighbour is no 
greater than half the expected separation value. This is 
to eliminate the contribution of the few points which are 
not part of the coherent flock (which does occur for the 
larger values of rj used in this study). Note our results are 
not sensitive to the choice of 7 (within a sensible range), 
indeed results are virtually identical if 7 G [0.25,0.75]. 
The a-shape is subsequently constructed from the points 
Yi, using a standard package in MATLAB. The results 
of this process can be seen in Fig. |31 This is then used to 
compute the morphological properties of the flock using 
Minkowski functionals, which we discuss next. 

Minkowski functionals measure the topological and ge¬ 
ometrical features of a given shape or structure. Had- 
wiger’s theorem [ 2 ^ states that in a d-dimensional space, 
there exist d -I- 1 quantities which completely describe 
the morphology of a given structure. In 2D, the three 
Minkowski functionals of a closed contour are its enclosed 
area S, perimeter P, and the Euler characteristic x [^ - 


but for the purposes of this study it is most convenient 
to analyse the dimensionless filamentarity 

F = {P^-47rS)/{P^ + ATrS) (7) 

By definition, F G [0,1], with F = 0 for a circle and 
F = 1 for a (not necessarily straight) line segment. Both 
S and P can be efficiently computed from the a-shape, 
and so we are able to compute the evolution of F in 
time, and report an average value for the hlamentarity, 
(F), as we did for Besag’s functions. Again we refrain 
from reporting results here, and present the bulk of re¬ 
sults once we have discussed a modification of the model. 
Note that the boundary dSa of the a-shape is a subset 
of the Delaunay triangulation of S, hence it is plausi¬ 
ble that we over estimate P and hence F due to small 
triangular artefacts as can be seen in Fig. [3] Again we 
emphasise that we are not attempting a detailed analysis 
of flock morphology, to compare with natural flocks. But 
to quantify the change in morphology when an external 
flow acts on the particles using a consistent methodology 
across all simulations. 


III. SELF-PROPELLED PARTICLES IN THE 
TAYLOR GREEN FLOW 

We now turn our attention to Vicsek flocks in the pres¬ 
ence of external noise, due to fluid motion. As discussed 
in section H] we forgo the computational expense, and 
complexity of a fully turbulent flow. Instead we turn to 
a steady vortical flow, the Taylor Green (TG) vortex 
defined as 

v/(x) = {uf,vf) = Vf{sm{x)cos{y),-cos{x)sin{y)), 

( 8 ) 

where x = (x, y). The vorticity field is given by 

a; = V X V/ = 2Vf sin(a;) sin(?/), (9) 

and the flow is incompressible (V-v/ = 0). V/ is a scaling 
parameter which can be adjusted to modify the relative 
intrinsic particle speed to that of the background flow. 
Figured! displays the vorticity field of the TG flow with 
the corresponding velocity field. 

A perhaps overly simplistic approach to modelling the 
system, would be to simply allow the particles to be ad- 
vected by the flow, such that the equation of motion for 
the particles is modified to become 

dx ■ 

= (V cos(6»i) -k Uf{xi),V sin( 6 >i) -k u/(xi)). (10) 

However, it is not clear to us how Eq. @ can still re¬ 
main valid. We would assume that ‘particles’ rather than 
orienting themselves to match the direction the other 
particles are oriented towards, they would instead orient 
themselves to the direction of motion of nearby particles. 
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FIG. 5: (Color online) Snapshots of the Vicsek flock with external fluid forcing, in the statistically steady-state (left rj = 
0.1, R = 0.2, right r) = 0.1, R = 0.8) with Vf = 0.5. A pseudocolour plot (made slightly translucent to improve visualisation of 
the particles) displays the flow’s vorticity uj, Eq. ([9ll. Arrows indicate the particles direction of motion. Note the concentration 
of particles in regions where oj ~ 0. 
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FIG. 6: (Golor online) The time averaged Besag’s function 
(L(r)), Eq. ([5} plotted as a function of spatial scale r for 
simulations with Vf = 0.5. L{r) > 0 indicates clustering of 
particles at a given scale, L{r) < 0 indicates segregation. 


Such an approach was also taken in m- Hence we arrive 
at a modified form for the equation used to update d, 

=atan2|^ ^ ^ cosdJ-fTyCz, 

\ |xi-Xj|<iJ \xi-Xj\<R ) 

where = atan 2 (Hsin( 0 j)-|-r;/(xj)), H cos( 0 j)-|-u/(xj)). 
We retain the intrinsic noise, to model the fact that it 
is unlikely real animals will perfectly align themselves 
with neighbours within the critical radius. As stated, we 
believe the justification of noise uncorrelated in space and 
time is reasonable to model intrinsic noise. Absent from 
this modified model is the rotation of particles due to 


the non-zero vorticity of the flow. The effect of particle 
rotation would make an interesting starting point to a 
future study, but we refrain from its inclusion in this 
work as it would require a substantial modification of 
the standard Vicsek model. 

Particles are evolved according to Eq. (HOD, with the 
particles direction determined using Eq. (EH)- Domain 
size and particle number remain unchanged. Q is drawn 
randomly at each timestep and we take At = 0.1. Unlike 
for the standard model, we use a 3’’'^-order Runge-Kutta 
time-stepping scheme, to allow for the additional com¬ 
plexity of the flow, however 9 is only updated at the 
beginning of a timestep. Again simulations (77 € [0,0.5], 
R € [0,1]) are evolved for 200 time units, and we monitor 
to ensure statistics are only computed in a statisti¬ 
cally steady value. We consider three values of U/ = 0.1, 
0.2 and 0.5, with a fixed value U = 1 for the particle 
speed. Snapshots of the particles for simulations with 
77 = 0.1, R = 0.2, 0.8, Vf = 0.5 are displayed in Fig. [5] 
along with the vorticity field of the TG flow. It is clear 
from a visual inspection, particularly for large values of 
R, that the morphology of the flock is dramatically af¬ 
fected by the influence of the flow. We note also that 
particles tend to be found in areas of low vorticity, this 
is not seen (and would not be expected) if we consider 
tracer particles U = 0. We shall discuss this feature later 
in the text. 

In all simulations we find the statistical properties of 
the flocks are nearly independent of magnitude of the 
intrinsic noise, and so for all quantities reported below 
we perform additional averaging over 77 . Figure | 6 ] shows 
the averaged Besag’s functions (Eq. ([S])) for five differ¬ 
ent values of interaction radius R, with Vf = 0.5. Again 
a noticeable (although less pronounced when compared 
with Fig. migration of the peak of (L) to the right is 
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FIG. 7: (Color online) (left) The time averaged characteristic size of the flock (D) (Eq. (HJ), plotted as a function of the radius 
of interaction R, for varying flow velocities Vf. (right) The corresponding plot of the variation of the flock hlamentarity (F) 
(Eq. ©). Note a profound change is observed in the flock morphology {R > 0.6) when Vf > 0.2. 



FIG. 8: (Color online) The probability density function 
(PDF) of (j), the angle between the flow direction and the 
particles intrinsic velocity, for the simulation with Vf — 0.2, 
rj = 0.1, R = 0.6 . The inset shows the PDF plotted with 
a linear scale. The main figure shows a logscale plot, with 
a power law fit PDF((()) ~ 0“^'^ displayed as a red (dashed) 
line. 


observed. However due to the periodic boundary condi¬ 
tions a secondary peak where {L{r)) > 0 for r > tt can 
be seen for all values of R. This indicates that the flow 
is forcing the flocks to form an elongated shape, tightly 
clustered in one direction, but coherent over the scale of 
the box, as we see in Fig. [S] 

This is also visible if we inspect plots of {D) (Eq. dH])) 
against R for the values of Vf considered. Fig. [7] (left). 
A clear reduction in the typical distance of clustering 
{D) is seen (at least for larger values of R) as the flow 
velocity increases. The effect of the flow on the flocks 
morphology, namely {F) (Eq. d?])), can be seen in Fig. [7] 
(right). At low flow speeds {Vf = 0, 0.1) increasing the 


radius of interaction leads to a dramatic decrease in the 
hlamentarity, as the flock tends to become close to cir¬ 
cular. However above some critical value of Vf, there is 
a dramatic transition, and even for large values of R the 
structures exhibit large hlamentarity. It is important to 
note that this is visible even when the how speed is only 
one-hfth of the particles intrinsic swimming speed. 

Finally we turn our attention to the relationship be¬ 
tween the particles intrinsic direction of motion, and the 
direction of the local how. We dehne 


Ipi = cos 


V ^|v/(x.0| / ’ 


( 12 ) 


from this we use kernel density estimation to esti¬ 
mate the probability density function (PDF) of (j). This 
is displayed in Fig. [51 We observe empirically that 
there is a power law relationship between the angle be¬ 
tween the how and the particles intrinsic velocity, i.e. 
PDF(0) ~ (/)“. Of course it is then natural to ask if this 
power law is universal or depends on rj, R and/or Vf. 
We address this question directly in Fig. [51 for the two 
larger how velocities Vf = 0.2, 0.5. It is clear that cj) is 
independent of the intrinsic noise ry, but depends on both 
the how speed Vf and the radius of interaction R. The 
behaviour is as one may naively expect with a much shal¬ 
lower power law relationship for large values of R, and 
smaller values of Vf, where one may expect more scat¬ 
ter between the how direction and the particles intrinsic 
direction. 


IV. DISCUSSION 

The main goal of this work was to understand the ef¬ 
fect of an external forcing how, made up of steady coun¬ 
terrotating vortices, on the standard Vicsek model. To 
do this we made a slight modihcation to the model and 
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FIG. 9: (Color online) The power law scaling, a, of the probability density function (PDF) of ifi (PDF(())) = <()“) where ifi is the 
angle between the flow direction and the particles intrinsic velocity, (left) Vf = 0.2, (right) Vf = 0.5, plotted for varying radii 
of interaction R. The slope is clearly dependent on both R and Vf, but appears independent of the intrinsic noise rj. 
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FIG. 10: (Color online) Snapshots of the Vicsek flock with 
external fluid forcing given by Eq. with rj — 0.1, R = 0.8 
and Vf — 0.5. A pseudocolour plot (made slightly translucent 
to improve visualisation of the particles) displays the flow’s 
vorticity uj. Arrows indicate the particles direction of motion. 


then compared statistically properties of the flocks, pay¬ 
ing particular attention to morphological properties. One 
of the striking features of the flocks with an external flow 
Vf > 0.2 is that there is a marked increase in the fila- 
mentarity of the flock, when compared to the zero flow 
case. As is clear in Fig.[S]this is due to the fact that par¬ 
ticles appear confined to areas of low vorticity. The phe¬ 
nomenon of spatial clustering of independent inertial par¬ 
ticles, i.e. subject to Stokes dragOn vortical and turbu¬ 
lent flows is a well explored field . Because the parti¬ 
cles posses inertia their flow is compressible, even though 
the background flow is typically incompressible. Hence 
large fluctuations of particle concentration are found, and 
particles tend to collect in elongated sheets (in 3D) on 


the edges of vortical structures. We see exactly the same 
phenomena here, not due to a Stokes drag, but because 
the intrinsic motion of the particles means they do not re¬ 
spond instantaneously to the flow. Hence they are ‘flung’ 
out of vortical regions and are concentrated in thin rib¬ 
bons outside of the vortical cells. This is seen to have a 
profound effect on the morphology of the flock. 

In order to investigate the universality of the results 
presented here we perform one further numerical simu¬ 
lation, with a random, unstructured flow. This is gener¬ 
ated from the summation of 3 random Fourier modes in 
the spirit of the KS model used in [HI- We take a 
velocity field defined as 

3 

^/(^) ~ ^ ^ (-^771 ^ kjTj COS 077 J Bttt, X kjTj, sin 0777 ), 

771 = 1 

(13) 

where 4>rn = k;,„ • x, are wave vectors, A,„ and 
are random unit vectors (orthogonal to k™) and v is a, 
scaling parameter chosen such that the root mean square 
(RMS) value of Vf is Vf. Taking km = w creates the 
random large scale (solenoidal) flow whose vorticity field 
is displayed in Fig. [Tni A simulation of the Vicsek flock 
in the presence of this random fluid flow is performed 
taking 77 = 0.1, R = 0.8 and Vf = 0.5; a snapshot of the 
system in the steady-state is also displayed in Fig. [TO] 
The elongated structure of the flock visible in the sim¬ 
ulations using the TG flow is also found in this random 
flow, which gives some indication of the robustness of the 
results presented here. 

Of course it is an open question as to the extent 
to which these results apply to real animal flocks and 
swarms. However it is both interesting and important 
to understand how structures in the flow can affect both 
the shape and fluctuations (which we have not consid¬ 
ered here) of the shape of animal flocks. As mentioned 
earlier in the text we have neglected how rotation of the 








particles affects their dynamics, which could provide an 
interesting extension to the model. Natural next steps 
are to understand the role that temporal fluctuations of 


the flow play, before turning attention to more realistic 
flows, such as those generated by direct numerical simu¬ 
lation of the Navier Stokes equation. 
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